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1.  INTRODUCTION 

To  minimize  clutter  problems  and  resolve  contacts  in 
ck  se  proximity,  harbor  surveillance  radars  should  use  high 
resolution  antennas.  The  influence  of  these  considerations 
on  the  choice  of  parameters  desired  for  the  two  radars  eventu¬ 
ally  installed  in  the  San  Francisco  Experimental  Vessel  Traf¬ 
fic  System  was  examined  during  the  early  development  activity. 
Although  the  antenna  specifications  projected  at  that  time  dif¬ 
fered  somewhat  from  those  finally  used,  the  Coast  Guard  has 
requested  that  the  basis  for  the  original  specifications  be  docu¬ 
mented;  that  is  the  intent  of  this  report. 

It  will  be  seen  that  clutter  considerations  dominated 
the  design  of  the  antenna  for  the  Point  Bonita  site  while  reso¬ 
lution  was  the  essential  constraint  on  the  Yerba  Buena  Island 
antenna. 


0 

11 

e 


T 


THC  JOHNS  HOPKINS  UNIVtKSmr 

APPLIED  PHYSICS  LABORATORY 
Sh.w*  $Mmt0  mastviamc 


2.  EA  JAR  RESOLUTION 

Three  sections  of  the  Yerba  Buena  Island  (YBI)  sur¬ 
veillance  area  place  restrictions  on  resolution:  the  entrance 
to  the  Oakland  Inner  Harbor  and  two  channels  near  Richmond 
(Ref.  1). 

In  the  former  case  (Fig.  1,  extending  radially  at  about 
4  o'clock),  the  junction  of  ♦•'a  Inner  Harbor  Entrance  Channel 
and  the  Inner  Harbor  Reach  narrows  to  500  feet  at  a  point 
1.  9  nmi  from  the  YBI  antenna.  Since  this  channel  is  approxi¬ 
mately  parallel  .o  the  radar  line  of  sight,  cross -channel  dis¬ 
crimination  must  be  accomplished  via  antenna  azimuth  reso¬ 
lution.  To  split  this  width  into  quarters,  the  required  beam- 
width  is: 

eOJ  tan  4-  1.9  *  6076 
=  0.620°. 

At  the  start  of  the  bend  further  into  the  reach,  the  channel 
width  is  still  only  600  feet;  however,  the  distance  from  the 
radar  has  increased  to  3.25  nmi,  and  the  resolution  require¬ 
ment  is  more  stringent: 

_  .  -1  600 
dOB  ”  an  4  •  3.25  •  6076 

=  0.435°. 

In  the  Richmond  area  (Fig.  2),  there  are  two  channels 
of  interest.  The  Point  Potrero  Reach  is  500  to  600  feet  wide, 
but  it  is  essentially  perpendicular  to  the  YBI  radial,  and 
cross-channel  discrimination  can  be  achieved  through  use  of 
a  narrow  radar  pulse  width.  The  northern  end  of  the  South¬ 
ampton  Shoal  Channtl  (11  o'clock  radial)  is  600  feet  wide, 

7.  2  nmi  from  the  YBI  site,  and  parallel  to  the  radar  radial. 

It  requires  an  o.zimuth  beamwidth  of; 

6  ,  tan"  1  _ 622 - 

SH  4  •  7.2  -  6076 

=  0.  196'. 
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Fig.  2  RICHMOND  AREA 
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Based  only  on  these  considerations,  the  Southamp¬ 
ton  constraint  would  define  the  azimuth  beamwidth  require¬ 
ment  for  the  YEI  radar.  However,  the  Coast  Guard  indi¬ 
cated  that  the  Oakland  Inner  Harbor  Reach  should  be  the 
basis  for  setting  the  azimuth  specification  on  the  grounds 
that  another  radar  is  likely  to  be  placed  at  Richmond  if  the 
system  is  ever  expanded.  To  resolve  the  Southampton 
Shoal  Channel,  the  Automatic  Detection  and  Tracking  por¬ 
tions  of  the  San  Francisco  system  use  data  received  on 
successive  dwells  to  determine  the  centroid  of  a  contact, 
thus  achieving  discrimination  between  contacts  separated 
by  less  than  the  bear  lwidth. 
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3.  SEA  CLUTTER 

Since  the  Point  Bonita  radar  looks  primarily  seaward 
(Fig.  3),  resolution  is  significant  only  in  separating  targets. 
The  2000  foot  wide  Main  Ship  Channel  is  more  easily  resolved 
than  the  Oakland  Inner  Harbor  Reach,  even  though  its  distant 
end  is  5.  8  nmi  from  the  Point  Bonita  site.  Sea  clutter  dom¬ 
inates  the  design  of  this  antenna.  A  pulse  radar  viewing  sea 
clutter  at  low  grazing  angles  has  a  maximum  range  given 
(Ref.  2)  by: 

_  2  o  cos  <t> 

max  (s/c)  a  ?  ct 
o 

where: 

a  =  target  cross  section, 

;  !  <fi  =  grazing  angle, 

s/c  =  minimum  signal-to- clutter  ratio, 
o0  =  sea  backscatter  coefficient, 
e  -  antenna  azimuth  beamwidth, 
c  =  velocity  of  propagation,  and 
,  T  =  radar  pulse  width. 

This  equation  can  be  employed  to  estimate  the  detectable 
•  ;  contact  cross  section  as  a  function  of  range  and  sea  state. 

(  |  However,  there  are  a  number  of  assumptions  which  influence 

the  result,  and  they  warrant  discussion, 
r 

|  The  first  problem  is  to  obtain  usable  values  from  the 

backscatter  coefficient.  We  have  employed  those  of  Ref.  3 
j  >  to  generate  Figs.  4,  5,  and  6.  The  tables  from  which  the 

[  !  data  points  were  obtained  were  compiled  using  experimental 

data  from  numerous  sources.  Assumptions  and  constraints 
are  summarized  in  the  reference.  For  our  purposes,  the 
'  ,■  most  serious  of  these  appears  to  be  that  the  data  was  based 

on  pulse  lengths  in  the  0.  5  to  10  us  region  with  echoes  hav¬ 
ing  approximately  Rayleigh  distributions.  But  the  VTS  x’adars 
will  have  50  ns  transmitters  so  that  clutter  distribution  can 
hp  pxnpcted  to  deoart  from  that  model. 

-  '  ir  -  i 

Despite  this  limitation,  the  data  has  been  employed 
in  the  hope  that  it  offers  some  guidance  in  the  selection  of  at 
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Fig.  4.  X-BAND  SEA  BACKSCATTER  COEFFICIENT,  HORIZONTAL  POLARIZATION 
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Fig.  5.  X-BAND  SEA  BACKSCATTER  COEFFICIENT,  VERTICAL  POLARIZATION 


Fig.  6. 
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least  antenna  polarization.  For  example,  consider  Fig.  7. 
Antenna  beamwidth  was  set  as  narrow  as  seemed  practical 
(i.  e. ,  0.  25°)  on  the  grounds  that  it  is  comparable  to  that  em¬ 
ployed  in  some  of  the  European  harbor  radars  and  leads  to  a 
not  unreasonable  aperture  on  the  order  of  30  feet.  Assuming 
a  16  dB  signal- to- clutter  ratio  is  required  to  prevent  false 
alarm  overload  of  the  radar  computer,  the  detectable  contact 
cross  sections  are  as  shown  for  the  two  polarizations.  Al¬ 
though  not  plotted,  horizontal  polarization  offered  an  even 
further  advantage  over  vertical  at  sea  state  1.  For  example, 
at  3  nmi,  the  detectable  cross  sections  were  0.029  and  0.  63 
m-  respectively.  Thus,  it  appears  that  horizontal  polariza¬ 
tion  might  be  the  better  choice. 

However,  the  literature  also  offers  evidence  that,  as 
pulse  lengths  decrease,  the  departure  from  Rayleigh  is  greater 
for  horizontal  polarization  than  for  vertical  (Refs.  2  through 
5).  This  might  prove  significant  enough  to  offset  an  apparent 
horizontal  polarization  advantage  for  the  following  reasons. 

The  requirement  to  minimize  radar  computer  satura¬ 
tion  sets  a  value  for  false  alarm  rate.  In  this  application,  it 
is  on  the  order  of  10"^.  Assume  we  require  a  0.  5  probability 
of  detection  and  are  unable  to  obtain  the  benefits  of  pulse-to- 
pulse  integration  (e.  g.,  because  of  the  high  resolution  and 
long  correlation  time).  If  vertical  polarization  yields  Ray¬ 
leigh  statistics,  the  required  s/c  ratio  is  10  dB  for  n  =  1. 

If  horizontal  polarization  results  in  log-normal  returns,  then 
to  achieve  the  same  false  aiarm  rate,  the  s/c  value  must  be 
much  greater  (e.g. ,  26  dB  for  sea  state  3).  Given  such  as¬ 
sumptions,  the  resulting  difference  in  detectable  contact 
cross  section  is  shown  in  Fig.  8. 

One  might  legitimately  challenge  the  specific  numeri¬ 
cal  difference  noted  in  the  figure  on  grounds  that: 

1.  Backscatter  coefficients  were  not  defined 

at  short  pulses  and  are  thus  suspect  to  some 
extent,  and 

2.  Vertical  polarization  has  some  log-normal 
tendencies  rather  than  being  pure  Rayleigh 
as  assumed  in  arriving  at  that  result. 
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NOTES: 

SEA  STATE  3 

FALSE  ALARM  RATE  -  10“6 
PULSE  WIDTH  -  SO  NS 
AZIMUTH  BEAMWIDTH  -  0.25* 
ANTENNA  HEIGHT  -  290  FEET 
PROBABILITY  OF  DETECTION  -  0.5 


HORIZONTAL  POLARIZATION,  LOG-NORMAL 
VERTICAL  POLARIZATION,  RAYLEIGH 


RANGE  (nmi) 


Fig,  8.  IMPACT  OF  ANTENNA  POLARIZATION  ASSUMING 
DIFFERENT  CLUTTER  DISTRIBUTIONS 
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However,  the  literature  does  offer  support  for  the  difference 
in  distributions.  At  higher  sea  states,  the  advantage  should 
be  further  in  favor  of  vertical;  with  calmer  waters,  the  ad¬ 
vantage  would  decline  or  ever  vanish.  However,  at  the  lower 
sea  states,  small  contact  detection  is  not  as  great  a  problem. 

Thus  we  conclude  that  the  Point  Bonita  antenna  should 
use  vertical  polarization.  For  the  calmer  conditions  en¬ 
countered  within  the  harbor,  horizontal  polarization  is  per¬ 
haps  the  better  choice,  since  the  calmer  waters  (in  which 
horizontal  polarization  has  an  advantage)  will  occur  a  grea¬ 
ter  portion  of  the  time  than  it  will  for  the  Point  Boniia  sur¬ 
veillance  area. 

If  the  absolute  values  in  either  Fig.  7  or  8  are  accu¬ 
rate,  one  could  argue  that  the  narrow  azimuth  beamwidth  con¬ 
stitutes  an  overkill  and  a  smaller  aperture  antenna  could  be 
deployed  for  Point  Bonita  at  less  expense.  But  the  many  un¬ 
certainties  associated  with  generation  of  the  curves  argues 
for  the  more  prudent  approach;  an  azimuth  beamwidth  on 
the  order  of  0.  25°.  The  YBI  azimuth  beamwidth  was  dictated 
by  resolution  considerations  as  0. 44°. 

In  light  of  experience  at  the  latter  site  with  rip  tides, 
vertical  (rather  than  horizontal)  polarization  would  probably 
have  been  a  better  choice,  although  we  have  no  clear  cut 
means  of  demonstrating  that  vertical  would  yield  the  same 
advantages  it  offers  against  conventional  sea  clutter. 
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4.  PRECIPITATION  CLUTTER 

Elevation  beamwidths  for  both  radars  are  determined 
by  requirements  for  performance  in  a  precipitation  environ¬ 
ment  and  surveillance  coverage  close  to  the  radar  sites.  \Ve 
will  consider  the  close-in  coverage  first.  Coverage  is  re¬ 
quired  to  within  about  1200  feet  of  both  the  PB  and  YBI  radars. 
To  determine  the  vertical  beamwidth,  we  first  assume  a 
Gaussian- shaped  antenna  pattern  (Ref.  6)  given  as: 

2 

f (a)  =  exp(-2.  88  -j-)  , 

where  4  =  the  3  dB  vertical  beamwidth 

a  =  the  angle  off  axis. 

The  bore  sight  of  the  antenna  is  depressed  such  that 
the  beam  illuminates  a  point  at  maximum  range  on  a  4/3 
earth.  From  Ref.  3,  this  is: 

q  .  -1/h  _R  . 

@  =  sm  <g  +  £->. 

e 

where: 

h  =  antenna  height 
R  =  slant  range,  and 

r  =  4/3  radius  of  the  earth  (4587). 
e 

Signal  strength  is  inversely  proportional  to  the  fourth  power 
of  range.  Thus,  for  any  given  beamwidth,  antenna  height, 
and  maximum  range  (290  feet  and  15  nmi  for  PB,  400  feet 
and  8  nmi  for  YBI),  one  can  estimate,  as  a  function  of  range, 
the  signal  strength  relative  to  that  at  maximum  range.  This 
is  plotted  for  the  two  radar  siies  in  Figs.  9  and  10.  To 
achieve  the  same  signal  strength  at  1200  feet  that  is  avail¬ 
able  at  maximum  range,  the  curves  indicate  that  the  verti¬ 
cal  beamwidth  for  the  Point  Bonita  antenna  should  be  approxi¬ 
mately  8°.  To  intercept  the  surface  at  maximum  range,  the 
boresight  should  be  depressed  0.  27°.  For  Yerba  Buena  Is¬ 
land,  the  corresponding  angles  are  11°  and  0.  52°,  respectively. 

These  elevation  beamwidths  combine  with  the  azimuth 
beamwidth  and  pulse  width  to  define  a  volume  in  space  which 
influences  the  level  of  precipitation  clutter  return.  As  a 
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Fig.  10.  INFLUENCE  OF  VERTICAL  BEAMWIDTH  ON  MINIMUM  RANGE,  YERBA  BUENA  ISLAN.O 
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result,  the  detectable  contact  cross  section  is  given  approxi¬ 
mately  by  Ref.  3  as: 

2 

a  =  R  (tt/ 8)  0 $ c  x  (s/c) 

where : 


R  =  range, 

0  =  azimuth  beamwidth, 

4>  =  elevation  beamwidth 
c  =  propagation  velocity, 

T  =  pulse  width, 

s/c  =  minimum  signal- to- clutter  ratio,  and 
E  =  precipitation  backscatter  coefficient. 


Applying  this  equation  using  backscatter  coefficients  of 
Ref.  3  yields  the  performance  for  the  two  radars  as  shown 
in  Figs.  1 1  and  12. 


The  solid  curves  are  for  linearly  polarized  antennas 
and  they  indicate  that  either  moderate  or  heavy  rain  can 
restrict  the  detection  of  2  m2  targets.  To  overcome  this 
limitation,  it  will  be  necessary  to  include  provisions  for 
selecting  circular  polarization  on  both  antennas  at  the  opera 
tor's  option.  Alternately,  one  could  consider  decreasing 
the  vertical  beamwidth  and  sacrificing  close-in  coverage. 
But  the  full  13  dB  improvement  obtainable  with  circular 
polarization  cannot  be  matched  without  a  drastic  reduction 
in  the  vertical  beamwidth. 


Figures  13  through  16  offer  performance  projections 
for  relaxed  azimuth  beamwidths. 
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Fig.  13.  YCR8A  BUENA  ANTENNA  OPTIONS 
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Fig.  15.  PRECIPITATION  CLUTTER  AT  POINT  BONITA 
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Fig.  16.  PRECIPITATION  CLUTTER  AT  YERBA  BUENA 
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5.  OTHER  CONSIDERATIONS 

Another  aspect  to  be  considered  in  this  antenna  de¬ 
sign  is  the  effect  of  lobing  on  the  detection  of  surface  con¬ 
tacts.  The  elevation  angle  of  the  first  lobe  is  approximately 
X/(4ha),  where  X  is  the  wavelength  and  ha  is  the  antenna 
height  (Ref.  5).  For  the  Point  Bonita  antenna  and  X-band 
operation,  this  is  (0.031  *  3.28l)/(4  •  290)  =  8.  5  •  10‘5 
radians  (z  0. 005°).  There  is  a  distance  at  which  the  free 
space  range  dependence  shifts  to  an  R"8  dependence  given 
approximately  as  (4haht)/X  ,  where  is  the  contact  height 
(Ref.  3).  Again,  for  Point  Bonita  and  X-band,  assuming 
a  3  foot  contact  height,  this  range  is 

(4  •  290  •  3)/(0.  031  •  3.28)  =  34  200  feet, 

or  5.  6  nmi.  For  this  same  antenna  and  S-band  operation, 
the  angle  and  range  are  0.  005  •  10/3. 1  =  0.016°  and 
5-  6  •  3.  1/10  =  1.73  nmi,  respectively.  This  last  result 
tends  to  make  S-band  operation  less  attractive  than  X-band 
because  of  its  impact  on  power  requirements.  For  the 
YBI  radar,  with  its  400  foot  antenna  height,  the  X-band 
transition  range  is  5.  6  •  400/290  =  7.7  nmi. 

The  distance  in  nautical  miles  to  the  radar  horizon, 
assuming  a  4/3  earth,  is  1.23  ha,  where  ha  is  in  feet  (Ref. 
4).  For  X-band  operations,  the  horizon  is  21. 0  nmi  and 
24.  6  nmi  for  PTB  and  YBI  respectively. 

Finally,  consideration  was  also  given  to  use  of  an 
inverse  cosecant  squared  vertical  beam  pattern  as  a  means 
of  achieving  close-in  coverage  and  obtaining  a  constant 
received  signal  independent  of  range.  The  possibility  has 
not  been  entirely  discarded  but  since  a  loss  in  gain  would 
result  and  the  broadened  beam  would  adversely  affect  pre¬ 
cipitation  performance,  a  decision  has  been  deferred  pend¬ 
ing  determination  of  the  remaining  radar  parameters.* 

The  CPS  software  employed  in  the  calculation  is 
included  as  Appendix  A  to  this  report. 

*The  final  antenna  design  did,  in  fact,  approximate  the  co¬ 
secant  squared  vertical  pattern  through  the  use  of  a  second¬ 
ary  feed. 
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APPENDIX  B 
CPS  SOFTWARE 


The  attached  listing  is  a  PL/I  program  used  to 
expedite  some  of  the  calculations  associated  with  creating 
the  figures.  It  is  used  at  a  conversational  terminal. 

For  example,  calling  CLUT  lets  the  user  enter  a 
signal-to-clutter  ratio,  pulse  width,  azimuth  beamwidth, 
and  antenna  height.  The  program  then  computes  the  graz¬ 
ing  angle  for  several  preset  ranges.  The  user  then  uses 
Fig.  4  or  5  to  obtain  the  backscatter  coefficient  for  each 
angle  and  enters  that  data.  The  program  replies  with  the 
resulting  detectable  target  cross-section,  thus  providing 
a  point  for  curves  like  those  in  Figs.  7  and  8. 

The  routine  PREC  will  provide  data  for  a  figure 
like  Fig.  15  when  the  user  supplies  arguments  for  azimuth 
angle,  elevation  angle,  pulse  width,  and  signal-to-clutter 
ratios. 
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OQC:  PROCEDURE  (a); 

c*TRUNC( a/ 1000 ) ; 

IF  o."0  THEN  b-a; 

If  c-0  THEN  GO  TO  Aj  ELSE  b"H0D(e,c); 

A:  e-SIGN(a); 

d"10«»( ABS(b)/10)*10"»c; 

IF  a>0  THEN  RETURN  (d); 

IF  e"Q  THEN  RETURN  <1); 

IF  a<0  THEN  RETURN  <  1/d >8 
END  OBC; 

SICC:  PR0CE0URE  (f,g); 

X*g/( f *6076  )-f / ( 2*6587 ); 

PS"ATAN(X,SQRT(1-X*»2>); 

PSD"PS"(  180/J. 16159); 

PUT  LIST('R"',f,'nm  Psl •' , PSD. 'degree^ 1 ) j 
END  SIGC; 

CLUTs  PROCEOURE  ; 

PUT  L I  ST  C 'S/C  ratio  (db),  Tau  (nsec),  Az  Beam  ioe*i,  Ant  Ht  (ft)') 
GET  LI STCSCO, TAU, A20,HF ) ; 

SC-DBC(SCD); 

TAS"TAU* . IE-08 ; 

AZR-AZO"(J. 16159/180); 

PUT  LIST{'  ‘); 

00  R"l, 1.5, 2, 2. 5, 3, 6, 6, 8, 10; 

CALL  SIGC(R,HF); 

PUT  LISTCSIGO  for  that  PS  I  * ) ; 

GET  LIST(SIGO); 

SIGR-UBC(SIGO); 

SIGTR-R*1852*S IGR*SC»AZR*. 3E J9 • . 5*TAS/C0SD ( PS  I ); 

PUT  LIST(R,'nm  SI GMA-', SI6 TR, *sq  m'); 

PUT  LIST'*  '); 

END  CLUT; 

ELAN:  PROCEDURE  (h,rmax, rmln); 

AA-SQRT(h«"2+rmln""2); 

B-SQRT ( h"»2*rmax*"2 ) ; 

C"rmax- rmln; 
aa*ATAN0(h, rmax); 

Sc*C*SIND(aa)/AA; 
cc-ATAN0( Sc, SQRT< 1-Sc**2 ) ) ; 

PUT  LISTt'El  beam  •' ,cc, 'decrees  1 ) ; 

dd*ATANO(rmln,h); 

ee*90-dd-cc/2; 

PUT  LIST< 'Depression  angle  »',ee, 'degrees' ); 

END  ELAN; 

ANT:  PROCEDURE  ; 

PUT  LIST('Max  Ground  Range  (nm).  Antenna  Ht  (ft)'); 

GET  LIST(rmax,h); 
rmaf *607b»rmax; 

PUT  LIST('Deamwldth  (degrees)'); 

GET  LIST(THET); 

RSMA"SQRT(h»»2*rmaf »*2 ) ; 

B£l"h/RSKA*RSMA/(2*6587«6076); 

B£TA"ATAND(BE 1, SQRT( 1-BE1*«2)); 

PUT  LI  ST ('BETA»', BETA); 

DO  ALF-0, 1,1. 5, 2, 2. 5, 3, 6, 6, I, 10, 15, 20; 

GAM"ALF*BETA; 

aq«l; 

bq"-2*6587*607E*SI NO(GAM); 
eq"2*h«6537»6076; 

RSL".5»(-bq-SQRT(bq»»2-6»aq«cq) ); 

GR"SQRT(RSL»*2-h»h ) ; 

IF  ALF"0  THEN  SSHR"1/RSL»«6; 

FG"EXP( -2.88«(ALF*»2/THET**2)); 
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136. 

SIGS-1/SSMR»FG»FG*( 1/RSL--6 ); 

131. 

POT  LISTL  *ALF-', ALF,  'GR-'.OP, 'SIC-*, LlOLsIGS), *db’ ); 

160. 

IF  SIGS<.9  THEN  GO  TO  JF; 

162 . 

ENO  ; 

166. 

JPt 

PUT  LISTC  * } i 

166. 

ENO  ANT; 

161. 

LET  L10(c>«10»L0G<c>/2. 302511) 

150. 

PREC: 

PROCEDURE  (az.al, pw.ae); 

152. 

PUT  LI ST( 'Output!  art  for  1,6,610  run'); 

156. 

PUT  Li  STL 'Fog,  drlz,  1«ht  mod,  hvy  (-,  .2S,1,6, 16  mm/hr)*>; 

156. 

00  ES— 100.-82,-72,-62,-53; 

156. 

PUT  LISTC  '); 

160. 

DO  R-l,6,10j 

162. 

RM-R-16S2; 

166. 

PI-3. 16159J 

166. 

S ICL»RH«RM*PI ’.25-az-al •, 3E09-1. 765 •, 000 1765* , 5*pw». lE-08«DBC(sc)*DIC(E5) 

161. 

PUT  LI  STL ' Clutt*r-',ES); 

170. 

PUT  LISTCR-'.R.'nm,  SIG-',SICL, 'tq  mtr*'); 

172. 

ENO  PREC; 

I 

* 


i 


